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e Adaption of coronaviruses OC43 and HKU1 to humans 
involved loss of HE lectin function 


e OC43 HE receptor binding site was lost via progressive 
accumulation of mutations 


e Loss of HE receptor binding alters sialate-9-O-acetylesterase 


receptor destroying activity 


e Balance of receptor binding and receptor destruction 
contributes to host tropism 
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In Brief 


Coronavirus OC43 entered the human 
population relatively recently. Bakkers 
et al. report that as an adaptation to 
replication in human airways, the OC43 
hemagglutinin-esterase lost its receptor- 
binding function. Consequently, virion- 
associated receptor-destroying activity 
toward clustered sialoglycan-based 
receptor determinants was reduced. 
Suggestive of convergent evolution, 
human respiratory coronavirus HKU1 
underwent similar changes. 
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SUMMARY 


Human betat-coronavirus (B1CoV) OC43 emerged 
relatively recently through a single zoonotic introduc- 
tion. Like related animal B1CoVs, OC43 uses 9-O-acet- 
ylated sialic acid as receptor determinant. B1CoV 
receptor binding is typically controlled by attach- 
ment/fusion spike protein S and receptor-binding/re- 
ceptor-destroying hemagglutinin-esterase protein 
HE. We show that following OC43’s introduction into 
humans, HE-mediated receptor binding was selected 
against and ultimately lost through progressive accu- 
mulation of mutations in the HE lectin domain. Conse- 
quently, virion-associated receptor-destroying activ- 
ity toward multivalent glycoconjugates was reduced 
and altered such that some clustered receptor popula- 
tions are no longer cleaved. Loss of HE lectin function 
was also observed for another respiratory human co- 
ronavirus, HKU1. This thus appears to be an adapta- 
tion to the sialoglycome of the human respiratory tract 
and for replication in human airways. The findings sug- 
gest that the dynamics of virion-glycan interactions 
contribute to host tropism. Our observations are rele- 
vant also to other human respiratory viruses of Zoo- 
notic origin, particularly influenza A virus. 


INTRODUCTION 


Coronaviruses (CoVs), long considered of high veterinary impact 
exclusively, are now generally recognized as zoonotic threats of 
pandemic potential in consequence of the 2002/2003 outbreak 
of severe acute respiratory syndrome (SARS) and the emer- 
gence of Middle East respiratory syndrome (MERS) in 2012 (de 
Wit et al., 2016). SARS-CoV was contained within 3 months after 
its discovery, while MERS-CoV causes a classical zoonotic 
infection with limited human-to-human spread and, as yet, is 
incapable of sustained community transmission (Reusken 
et al., 2016; Zumla et al., 2015). However, four other respiratory 
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coronaviruses—alphacoronaviruses NL63 and 229E and lineage 
A betacoronaviruses OC43 and HKU1—successfully breached 
the species barrier and are currently maintained in the human 
population worldwide through continuous circulation (Su et al., 
2016; de Groot et al., 2011). Conceivably, the study of these 
genuine human coronaviruses (HCoVs) may yield clues to what 
is required for viral adaptation to the human host and thereby in- 
crease our understanding of the probabilities and risks of coro- 
navirus cross-species transmission. 

OC43 and HKU1, while generally associated with benign com- 
mon colds in healthy immunocompetent individuals, may cause 
significant morbidity and even mortality in the frail (Morfopoulou 
et al., 2016; Woo et al., 2005a). OC43, the best-studied HCoV, 
apparently arose relatively recently, 120 to 70 years ago, with 
the most recent common ancestor of all extant OC43 variants 
dating to the 1950s (Lau et al., 2011; Vijgen et al., 2005, 2006). 
OC43 groups in the species Betacoronavirus-1 (B1CoV), 
together with highly related viruses from ruminants (bovine coro- 
navirus, BCoV), swine (porcine hemagglutinating encephalomy- 
elitis virus, PHEV), equines (equine coronavirus, ECoV), leporids 
(rabbit coronavirus HKU14, RbCoV), and canines (canine respi- 
ratory coronavirus, CRCoV) (de Groot et al., 2011; Figure 1A). 
The extraordinary radiation of B1CoVs might be explained from 
their receptor usage, as they attach to 9-O-acetylated sialic 
acids (9-O-Ac-Sias) (Matrosovich et al., 2015), i.e., glycan com- 
ponents common in mammals and birds (Traving and Schauer, 
1998). Paradoxically, however, most B1CoVs, including OC43, 
have very narrow host ranges. They form distinct monophyletic 
clades congruent with host selectivity, and phylogenetic evi- 
dence strongly argues against recurrent inter-species transmis- 
sions (Lau et al., 2011; Vijgen et al., 2006; see also Figure 1A). 
With the emergence of OC43 seemingly sparked by a singular 
one-time zoonotic event, the founder virus likely possessed 
unique traits already to allow efficient infection of, and transmis- 
sion among, humans. In turn, these traits and subsequent 
adaptations to the new niche must have closed the door on rein- 
troduction of the virus into animals. Despite the close genetic 
relationship between OC43 and BCoV-Mebus, particularly 
(96.6% identity across their entire genomes), their high-preva- 
lence endemicity/enzooticity, and the scale and frequency of 
interactions between their host species, there is no evidence of 
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Figure 1. Loss of HE-Mediated Receptor 
Binding in Human Betacoronaviruses 

(A) Evolutionary relationships among lineage A be- 
tacoronaviruses. Neighbor-joining tree based on 
BCoV lineage A ORF1b sequences in the NCBI 
database (n = 206), with 100% bootstrap support 
for all major branches. Evolutionary distances were 
computed using the Maximum Composite Likeli- 
hood method in MEGA6 (Tamura et al., 2013). The 
positions of human coronaviruses OC43 and 
HKU1 are highlighted relative to those of animal 
lineage A betacoronaviruses, including the various 
B1CoV subspecies and classical mouse hepatitis 
virus-type murine coronavirus (MuCoV). ChRCoV- 
HKU24, Chinese rat coronavirus HKU24 (Lau et al., 
2015); LAMV, Longquan Aa mouse coronavirus 
(Wang et al., 2015); LRLV, Longquan RI rat coro- 
navirus (Wang et al., 2015). The position of CRCoV, 
a recent split-off of BCoV, is indicated by a red dot. 
(See also Figure S1.) 

(B) HE°-Fc lectins (2-fold serial dilutions, starting at 
50 ng/L) were compared by sp-LBA for relative 
binding to BSM (with 50 ng/L BCoV-HE°-Fc set 
at 100%). 

(C) Conventional HAA with 2-fold serial dilutions of 
HE°-Fe fusion proteins of B1CoV members and 
HCoV-HKU1 (starting at 25 ng/well). Wells positive 
for hemagglutination are encircled. 

(D) High-sensitivity nanobead HAA. Non-com- 


PHEV Ov HKU @OOOOODBDBIOWOE!S plexed nanobeads (“beads only”) and nanobeads 
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OC43 spreading to cattle nor of BCoV (or any other B1CoV) 
efficiently spreading into humans (Figure 1A). 

B1CoV attachment involves two distinct types of surface pro- 
jections: large 20-nm “spikes” comprised of homotrimers of the 
class | fusion protein S and stubby 8-nm protrusions that are ho- 
modimeric assemblies of the hemagglutinin-esterase protein HE 
(de Groot et al., 2011). The S protein, common to all corona- 
viruses, mediates receptor binding and fusion of the viral and 
host membranes (Heald-Sargent and Gallagher, 2012). HEs are 
found in toroviruses and in influenza C and D viruses, but among 
coronaviruses, only in lineage A betacoronaviruses (de Groot, 
2006; de Groot et al., 2011; Hause et al., 2014; Matrosovich 
et al., 2015). HE monomers have a bimodular structure with a 
carbohydrate-binding (“lectin”) domain appended to an enzy- 
matically active sialate-O-acetylesterase (“esterase”) domain 
(Langereis et al., 2009; Rosenthal et al., 1998; Zeng et al., 
2008). Typically, in B1CoVs, both the spikes and HEs bind 
9-O-Ac-Sias, while the HE esterase domain promotes virus 
elution through receptor destruction. In consequence of the 
opposing activities of receptor binding and receptor destruction, 
B1CoV attachment to sialoglycans is dynamic and reversible. 
Thus, dead-end binding of virions to decoy receptors in oropha- 
ryngeal, respiratory, and gastrointestinal mucus may be pre- 
vented. Moreover, in infected host cells, HE-mediated receptor 
destruction is essential for efficient release of viral progeny (Des- 
forges et al., 2013). 

Here we present a comprehensive structure-function study of 
OC43 HE. We demonstrate that over decades after OC43’s intro- 


duction, its evolution was marked by a progressive loss of HE re- 
ceptor-binding activity through the accumulation of select muta- 
tions in the HE lectin domain. The effect of these mutations on 
the organization of the carbohydrate-binding site (CBS) and on 
receptor binding is explained from the crystal structure of the 
BCoV HE-receptor complex (Zeng et al., 2008) and visualized 
directly by the structure of a contemporary OC43 HE solved to 
2.45A resolution. Evidence is provided that loss of HE recep- 
tor-binding activity resulted in a reduction of virion-associated 
esterase activity toward multivalent clustered substrates. We 
propose that inactivation of the HE lectin domain altered the 
balance between virion binding and esterase-mediated virion 
elution, apparently as an adaptation to the sialoglycome of the 
human respiratory tract, and we speculate that this may be a 
contributing factor to host selectivity. This view is supported by 
our observation that the HE of HCoV-HKU1 also lost its recep- 
tor-binding properties. A mechanism is proposed in which 
accessibility of receptors to destruction depends on HE lectin 
function in relation to S-HE size differences and in which the bal- 
ance between attachment and catalysis-driven virion elution is a 
determinant of host tropism. 


RESULTS AND DISCUSSION 


Lectin Properties of B1CoV HEs 

A comprehensive set of B1CoV HEs, expressed as esterase- 
inactive Fc-fusion proteins (HE°-Fc), was compared for their 
Sia-binding properties. In accordance with previous findings 
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(Langereis et al., 2015), the HEs of most animal B1CoVs bound to 
bovine submaxillary mucin (BSM) in a 9-O-Ac-Sia-dependent 
fashion in solid-phase lectin-binding assays (sp-LBA; Figure 1B). 
Remarkably, however, those of PHEV strain VW572 and proto- 
typic OC43 strain USA/1967 (also known as ATCC-VR759; Mclin- 
tosh et al., 1967) showed no detectable binding. In hemaggluti- 
nation assays (HAA), more sensitive than sp-LBA, PHEV HE° 
tested positive, albeit weakly (Figure 1C). OC43 USA/1967 
HE®°, however, did not hemagglutinate. To augment HAA sensi- 
tivity through multivalency-driven high-avidity binding, we com- 
plexed HE°-Fc chimeras to protein A-coated nanobeads (Fig- 
ure 1D). For PHEV HE®, sensitivity was increased 250-fold 
as compared to the standard assay. Under these condi- 
tions, modest, but reproducible, hemagglutination by OC43 
USA/1967 HE° was detected (Figure 1D). Apparently, it has 
lost most, though not all, of its lectin function. 

Our observations prompted the question of whether loss of HE 
lectin activity is a strain-specific trait resulting from adaptation to 
in vitro propagation or a characteristic also shared by naturally 
occurring OC43 viruses. We therefore RT-PCR-amplified HE 
genes from more recent sputum-derived OC43 strains (respira- 
tory season 2005). The encoded proteins OC43/NL/A/2005 
and OC43/NL/B/2005 HE®, representative for the two major HE 
lineages in OC43 (designated “A” and “B,” Figure S1 and vide 
infra), did not show any detectable binding by sp-LBA or 
nanobead HAA (Figures 1B-1D). Apparently, loss of affinity for 
9-O-Ac-Sias as seen for OC43/USA/1967 HE is not an artifact. 
In fact, the data suggest that HEs of contemporary OC43 vari- 
ants have lost receptor-binding activity altogether. 


Lectin Function of OC43 HE Impeded by a Combination 
of Mutations 
The HE ectodomains of BCoV strain Mebus and OC43 
USA/1967, each 365 residues in length, differ at 18 positions 
only (Figures 2A and 2B). To identify the differences accounting 
for loss of lectin activity in OC43 USA/1967 HE, we introduced 
OC43-specific substitutions in BCoV-Mebus HE°-Fc, in sets 
and individually, and tested for a reduction in lectin activity by 
HAA (Figure S2A). Ina complementary approach, OC43 HE res- 
idues were systematically replaced by BCoV HE orthologs (Fig- 
ures S2B-S2D). The results show that the mutations in the 
esterase and MP domains of OC43 USA/1967 HE (Figures 2A- 
2C; Figure S2A) do not affect ligand binding. In fact, loss of 
9-O-Ac-Sia binding can be attributed to four out of eight substi- 
tutions in the lectin domain (T'"4N, R'”’P, E'”8Q, and F24L; Fig- 
ures S2A-S2D). Combined replacement of these residues in 
OC43 USA/1967 HE® by BCoV orthologs restored binding affinity 
to that of BCoV-Mebus HE (Figure 2C). Individual replacement of 
these residues in the context of a restored OC43 USA/1967 HE 
(“TREF”) showed that each mutation affects receptor binding 
to more (T'"4N, E'’8Q) or lesser (R'”’P, F?*’L) extent (Figure 2C). 
The impact of these mutations on 9-O-Ac-Sia binding can be 
understood from the crystal structure of the BCoV-Mebus 
HE-receptor complex (Zeng et al., 2008). The BCoV HE recep- 
tor-binding region is comprised of six surface loops, five grafted 
on the beta-sandwich core of the lectin domain and one 
emanating from the esterase domain (Figure S3A). The actual 
carbohydrate-binding site consists of a deep hydrophobic 
pocket P1 and a more shallow hydrophobic depression P2 
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that accommodate the methyl groups of the Sia-9-O- and 
—5-N-acetyl moieties, respectively (Figure 4E). The side chain 
of F*"" (a residue in the 812/813 B-hairpin) separates P1 from 
P2 and, in the HE-receptor complex, intercalates between the 
Sia-acetyl groups. Ligand binding, thus largely based on shape 
complementarity and hydrophobic interactions, is stabilized by 
extensive protein-sugar hydrogen bonding, particularly involving 
the B7-B10 B-hairpin, through main chain atoms of L?'? and 
N?'4 to the sialate-5-N-acyl amide and the Sia carboxylate, 
respectively, and through the side-chain hydroxyl of S?'° to the 
sialate-C8-hydroxyl and -C9-acyl oxygen. Y'®* in the 85-66 
loop (residues 176-185) is particularly important for ligand bind- 
ing as its side-chain hydroxyl group hydrogen bonds with the 
sialate-9-O-acetyl carbonyl, while its aromatic ring together 
with the side chains of F2"', L?°°, and L?°” walls the P1 pocket 
that is key to receptor recognition (Figure 4E). 

Residues of the «5-a6 loop (residues 111-118) do not directly 
contact the ligand. However, the T''“N substitution in OC43 
HE created a novel N-linked glycosylation site (TTS — NRS) (Fig- 
ure 2A). Aglycan attached to this site would extend into the CBS 
and potentially reduce receptor binding through steric hin- 
drance. Indeed, the introduction of this glycosylation site in 
BCoV HE? resulted in an almost complete loss of binding as 
detected by sp-LBA (Figure 2C) and a 250-fold loss of apparent 
binding affinity as measured by HAA (Figure 2D; Figure S2A). 
The same mutant protein, but now expressed in N-acetyl 
glucosamine transferase | (GnT1)-deficient HEK293S GnTI- cells 
instead of HEK293T cells, showed a less pronounced loss of 
binding affinity (30-fold in HAA; Figure 2D). The data show that 
glycosylation at N'"* indeed reduces the affinity of the HE lectin 
domain for 9-O-Ac-Sias and that this reduction is positively 
correlated with glycan size/complexity. 

In BCoV HE, E'”® fixates the B5-86 loop by engaging in double 
hydrogen bonds with the backbone nitrogen atoms of S'°° and 
A'®® in the B4-85 loop (Figure S3B). We offer that in OC43 HE, 
the replacement by Gln and consequential loss of inter-loop 
hydrogen bonding destabilizes the B5-86 loop and that its 
increased flexibility impairs the Sia-binding site, presumably by 
affecting the critical P1 pocket. The consequences of the other 
two substitutions in OC43 USA/1967 HE are less evident, 
although it can be envisaged that the R'’’P replacement might 
again affect folding and/or stability of the 65-66 loop and 
thereby indirectly the positioning of the critical Y'®* residue. 
The F*’L substitution is more difficult to explain as the affected 
residue locates to the core of the HE lectin domain, distal to the 
CBS, and the mutation apparently decreases ligand binding 
through indirect long-range effects. The combined data show 
for the early OC43 isolate USA/1967 that the reduced lectin ac- 
tivity of its HE resulted not from a single but from a combination 
of mutations. 


Progressive Loss of HE Lectin Function during OC43 
Evolution 

In the course of this study, many full genome sequences from 
OC43 field variants from the US, Western Europe, and China, 
well documented with respect to place and date of virus sam- 
pling, became available in GenBank. This wealth of information 
allowed us to study the recent evolution of the OC43 HE protein. 
In phylograms, in accordance with the presumptive date of 
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Figure 2. Loss of OC43 USA/1967 HE? Lectin Affinity Attributed to a Combination of Four Mutations 

(A) Alignment of BCoV, OC43 USA/1967, OC43 NL/A/2005, and OC43 NL/B/2005 HE with sequences color coded by domain (membrane-proximal domain, red; 
esterase domain, green; lectin domain, blue). Residues crucial for esterase activity (SGNDH) are annotated. Amino acid differences are marked in black. 
N'"* substituting for Thr is colored in cyan, and the resulting N-linked glycosylation site (NRS) marked with red dots. D?° is marked with a red asterisk. 

(B) Overall structure of BCoV HE (PDB: 3CL5) with mutations in OC43 USA/1967 HE indicated. Domain coloring as in (A). 

(C) Comparison of binding affinities of BCoV HE°, OC43 USA/1967 HE®, and derivatives by sp-LBA. BCoV, BCoV-Mebus HE®; OC43 TREF, OC43 USA/1967 HE® 
with N'"4, p1’7, Q178 and L?4” replaced by BCoV-Mebus HE orthologs; T'“N, R'”’P, E'”°Q, and F*4”L, TREF derivatives with indicated residues re-converted to 
the autologous OC43 orthologs. 

(D) Conventional HAA with BCoV-Mebus HE? T'"4N expressed in HEK293T or HEK293S GnTI° cells. (See also Figure S2.) 
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Figure 3. Complete Loss of HE Lectin Func- 
tion during OC43 Evolution due to Progres- 
sive Accumulation of Mutations 

(A) Close up of the BCoV-Mebus HE CBS in 
complex with a-Neu5,9Ac22Me (in sticks). Resi- 
dues comprising the metal-binding site (MBS) are 
also shown in sticks, the potassium ion is shown 
as a magenta sphere, and interactions between 
the metal ion and coordinating amino acid resi- 
dues as red dashed lines. 

(B) Disruption of the MBS in BCoV HE leads to loss 


@ BCoV 
® BCoV D220y 
A BCoV D220H 


P oP bos 
v © 8 se of lectin affinity. sp-LBA as in Figure 1B. 
= a HE®-Fe (ng/ul) (C) Disruption of the BCoV HE lectin MBS 
through D?2°Y substitution renders receptor 
c Nanobead HAA D Standard HAA as ar even by high-sensitivity 
nanobeat 7 
BCoV 900000000008 Beow-Mebts 200000000006 (D) The Y*?°D mutation partially restores lectin 
BCov DY GOOG SOs OO i ae affinity of OC43 NL/B/2005 HE® 
0043 NLU/B/2005 20 OOO aear 7 eee ears . . 
A Nisnobead WAR (E) D°*"H substitution in BCoV HE results in ther- 
HE®-Fe anones molability of receptor binding. Conventional 
BCoV-Mebus 909009000000 HAA before (4°C) and after (87°C) a temperature 
OC43 NL/B/2005 ©) = & shift up. 
e DEsS NHBIZN0E END 00000000. (F) Receptor binding of OC43 NL/A/2005 HE® and 
: _ Standard HAA Se derivatives as determined by sp-LBA. In OC43 
BCoV pacactctctet a 4+ 4 | aE HE®-Fe NL/A/2005 HE®, N‘"4, P17”, Q'7®, and L?4” were 
Bcov pH HOOOO00 ) replaced by BCoV-Mebus HE orthologs (TREF), in 
Heh + <6°4:4° 4°44. 4-4-4 4 0 aR combination with (1) H?2°D substitution (TREF, 
BCoV D2°H OOO IOOO H??°D), (2) repair of the 85-B6 loop by substituting 
es D'®yyy"® for IIT (TREF, —IIT), or (3) H?2°D and 
HE°-Fc repair of the B5-B6 loop (TREF, H?°D, —IIT). 
F 
100 
iS @ BCov ; 7 
S 75 wi charge of the metal ion. In BCoV HE’, 
£ * Tees He 3 the D?2°Y mutation resulted in loss of all 
= 50 3 TREF, -lIT z detectable binding to 9-O-Ac-Sia (Fig- 
= - A TREF, H22, IT |Z ures 3B and 3C). In turn, a Y?2°D back 
Ss mutation in the type B HE of OC43 


HE°-Fc (ng/pl) 


OC43 emergence, the HE of the USA/1967 strain is placed close 
to the root (Figure S1). In the decades thereafter, the HE proteins 
divided into three clades. Clade C, the least populous, is 
comprised of HEs from OC43 strains sampled in the early 
1990s with no recent representatives and may have been re- 
placed by clades A and B. Comparative sequence analysis in 
combination with structure-function analysis showed that both 
type Aand B HEs acquired several mutations in the lectin domain 
additional to those already present in OC43 USA/1967, which 
explains the loss of HE receptor binding in extant OC43 variants. 
Unexpectedly, in B-type HEs, but not in those of types A and C, 
glycosylation at N''4 was lost again due to a S''®A substitution 
(NRS—NRA). Further inspection revealed the replacement in 
B-type OC43 HEs of D®”° by Y. D?2° is part of a potassium ion- 
coordinating metal-binding site (MBS), a signature element of 
coronavirus HE lectin domains involved in the organization of 
the Sia-binding site (Figure 3A) (Zeng et al., 2008). The MBS is 
formed by main-chain carbonyl atoms of $72', E76°, and L?”, 
together with the side chains of Q?2?, S765, and D2°, the latter 
of which is particularly important as it balances the positive 


360 Cell Host & Microbe 27, 356-366, March 8, 2017 


NL/B/2005 restored Sia-binding activity 

to levels detectable by high-sensitivity 

nanobead HAA (Figure 3D). Again, the 

structure of the BCoV HE-receptor com- 

plex provides an explanation for the 
detrimental effect of the mutation. As L?6° and L?°’ are in a 
loop that is stabilized through the coordination of the potassium 
ion, disruption of the MBS might affect the folding of this loop, 
the positioning of the Leu residues, and thereby the structure 
of receptor-binding pocket P71. 

Remarkably, also in A-type OC43 HEs, the MBS was targeted 
by substitution of D2?°, now by His (Figure 2A). As evident from 
comparative sequence and phylogenetic analysis, the substitu- 
tions in A- and B-type HEs were independent events, and their 
selection must be ascribed to convergent evolution. Again, 
D??°H mutation in BCoV HE resulted in loss of detectable Sia 
binding as measured by sp-LBA (Figure 3B). Unexpectedly, 
however, in conventional HAA routinely performed at 4°C, 
wild-type and mutant BCoV HEs hemagglutinated to similar 
extent (Figure 3E). The effect of the mutation became evident af- 
ter a shift up to 37°C, upon which hemagglutination with wild- 
type BCoV HE remained stable while that of the mutant protein 
resolved (Figure 3E). Apparently, substitution of D?2° by His is 
less disruptive than by Tyr, but results in thermolability of HE re- 
ceptor binding, likely to cause reduced binding affinity under 


physiological conditions. Repair of the OC43-specific substitu- 
tions already present in the HE of the 1967 strain (Figure 2A), in 
conjunction with a H22°D back substitution, did not restore lectin 
activity of OC43 NL/A/2005 HE (Figure 3F), suggesting the pres- 
ence of at least one additional mutation to prevent receptor bind- 
ing. Comparative sequence analysis revealed a major change 
uniformly shared by A-type HEs again in the B5-B6 loop, with 
four adjacent residues 183-186, Asp-Tyr-Tyr-Tyr, replaced by 
lle-lle-Thr (Figure 2A), apparently as a result of a double frame- 
shift mutation (Figure S3C). Repair of this mutation, in combina- 
tion with restorative changes at the other five sites, was required 
for OC43 NL/A/2005 HE to regain receptor-binding activity (Fig- 
ure 3F). The impact of the 65-86 loop mutation is easily under- 
stood from the structure of the BCoV HE-receptor complex, as 
it involved loss of Y'®*, an essential residue for ligand binding. 


Crystal Structure of an OC43A-type HE 
To assess the consequences of the combined mutations in the 
type A HE lectin domain, we determined the crystal structure of 


Figure 4. Crystal Structure of OC43 


NL/A/2005 HE 

(A and B) Side-by-side cartoon representations of 
the overall crystal structures of BCoV-Mebus HE 
(A) and OC43 NL/A/2005 HE (B). HE monomers are 
colored gray or by domain (as in Figure 2A). 

(C) BCoV and OC43 HE have identical esterase 
catalytic sites. Overlay of BCoV (blue) and OC43 
NL/A/2005 (green) HE esterase domains. Cartoon 
representations with residues crucial for activity 
indicated as sticks. 

(D) Overlay of BCoV and OC43 NL/A/2005 HE 
lectin domains. F*"' indicated with sticks. 

(E) 9-O-Ac-Sia (with carbon atoms in cyan) binding 
in the BCoV-Mebus HE lectin CBS as observed in 
the crystal complex (PDB: 3CL5). Close up with 
contacting amino acid side chains shown in stick 
representation. Hydrogen bonds are shown as 
black dashed lines, and hydrophobic interactions 
with the Sia-9-O- and —5-N-methyl groups as thin 
gray lines. P1 and P2 indicate the pocket and hy- 
drophobic depression, which accommodate the 
methyl groups of the Sia-9-O- and —5-N-acetyls, 
respectively. 

(F) Close up of the inactivated CBS of OC43 
NL/A/2005 HE. Residues corresponding to those 
in (E) are in stick representation and colored by 
atom type. NAG, N-acetylglucosamine attached to 
N''*, is shown in stick representation. (See also 
Figure S3 and Table S1.) 


OC43 NL/A/2005 HE at 2.45A resolution 
(for crystallographic details, see Table 
S1; PDB: 5N11). In overall fold, OC43 
HE is highly similar to other coronavirus 
HEs (Figure 4). Differences in the 
esterase domain are limited to substitu- 
tion of five surface residues, none likely 
to impair esterase function or activity. 
OC43 HE is also identical to BCoV HE 
in most of its lectin domain, particularly 
in the beta-sandwich core, but the loops 
shaping the receptor-binding site have changed dramatically 
(Figure 4D; Figure S3D). The 811-612 and 85-86 loops are 
disordered, suggestive of extensive flexibility, in the case of 
the latter loop in accordance with loss of E'’®-mediated inter- 
loop hydrogen bonding. Those (parts of the) loops that can be 
visualized are reoriented and displaced with respect to their 
original position in BCoV HE, utterly destroying the CBS. For 
example the 67/610 B-hairpin, comprising F2"', L?'7, S?'S, and 
N?"4 that are key to protein-ligand interaction in BCoV HE 
(Zeng et al., 2008), has shifted by approximately 7A (Figure 4D). 
The P1 pocket, arguably the most critical element of the CBS, 
no longer exists as (1) F2"' is no longer at its original position 
and its side chain is rotated by 90°, (2) Y'®* was lost as a result 
of the frameshift mutation in the B5-B6 loop, and (8) the side 
chains of L266 and L?°’ are reoriented apparently due to the 
loss of the MBS (Figures 4E and 4F). Finally, attached to N‘''4 
is an N-acetylglucosamine, providing formal evidence that the 
newly introduced glycosylation site in OC43 HEs is functional. 
Although the remaining residues of the sugar chain could not 
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be visualized, the glycan would stretch across the remains of 
the CBS (Figure 4F). 


Loss of HE Receptor Binding in HCoV-HKU1: A Case of 
Convergent Evolution? 

OC43 and HKU1 occupy the same niche, prompting the ques- 
tion whether they have been subject to similar selective con- 
straints. HCoV-HKU1 is also a lineage A betacoronavirus but 
distantly enough related to the B1CoVs to be assigned to a sepa- 
rate species (de Groot et al., 2011) (Figure 1A). Like OC43, 
HCoV-HKU1 binds to 9-O-Ac-Sia receptor determinants via its 
spike (Huang et al., 2015), and it possesses an HE (Woo et al., 
2005b). The HE of the prototypical HCoV-HKU1 2005 Hong 
Kong genotype A strain (Woo et al., 2006) tested negative for 
9-O-Ac-Sia binding by sp-LBA, standard HAA, and nanobead 
HAA (Figures 1B-1D). Comparative sequence analysis, including 
genotypes A and B, of HCoV-HKU1 (Woo et al., 2006) shows that 
the HKU1 HE has undergone massive deletions, as a result of 
which most of the lectin domain was lost (Figure S4). 


Loss of HE Lectin Function Reduces Receptor- 
Destroying Enzyme Activity toward Multivalent 
Substrates 

In influenza C and D viruses, the hemagglutinin-esterase fusion 
protein is uniquely responsible for viral attachment and entry, 
and hence HE function is receptor binding first and foremost. 
The same holds for the HEs of MHV-type murine betacoronavi- 
ruses (MuCoV) that mediate virion binding to O-Ac-sialoglycan- 
based attachment factors (Langereis et al., 2010), while S binds 
to the proteinaceous entry receptor CEACAM-1 (Peng et al., 
2011). At the other end of the spectrum, in OC43 and HKU1, 
HE’s sole remaining function is receptor destruction, with virion 
attachment to 9-O-Ac-sialoglycans exclusively assigned to S. 
For the animal B1CoVs that retain a functional HE lectin domain, 
however, the role of HE is less obvious. On the one hand, HE 
might participate in attachment to increase virion binding avidity. 
On the other, its lectin domain may serve primarily to promote 
catalytic activity toward high-multivalency substrates by 
bringing the esterase in close proximity to and prolonged asso- 
ciation with clustered glycotopes, in analogy with the carbohy- 
drate-binding modules of cellular glycoside hydrolases (Boras- 
ton et al., 2004). In the latter case, the HE lectin CBS in the 
human lineage A BCoVs may have been lost to downregulate re- 
ceptor-destroying activity. To study this, we compared soluble 
recombinant HEs of BCoV, OC43, and HKU1 for their relative 
enzymatic activity toward the monovalent substrate p-nitrophe- 
nol acetate (pNPA) and the multivalent substrate BSM. The latter 
glycoconjugate carries hundreds of O-linked glycans clustered 
in such densities to give the protein a filamentous bottlebrush 
appearance typical for mucins (Zappone et al., 2015). Whereas 
all three HEs showed comparable specific activities when as- 
sayed with pNPA (Figure 5A), those of OC43 and HKU1 were 
more than 250-fold less active than BCoV HE in receptor- 
destruction assays with BSM (Figure 5B). Inactivation of the 
BCoV HE lectin CBS through a F2''A substitution reduced its 
esterase activity toward BSM to that of wild-type OC43 HE. 
Conversely, restoration of the OC43 HE lectin domain increased 
esterase activity to that of BCoV HE (see “TREF,” Figure 5B). 
Thus, BCoV and HCoV HEs indeed differ in their reactivity toward 
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clustered receptors in consequence of the respective conserva- 
tion or loss of lectin function. 


Consequences of Loss of Lectin Function for Virion- 
Mediated Receptor Destruction 

OC43 HE efficiently depletes endogenous pools of O-Ac-Sias 
within infected cells (Figure S5) as is critical for the release of viral 
progeny (Desforges et al., 2013). The loss of HE lectin function 
therefore is more likely to affect virus biology at another stage, 
namely that of virion (pre)attachment. Of note, the observations 
for recombinant soluble HEs cannot be extrapolated one-to- 
one to virion-associated HEs, because in virus particles, HE 
esterase activity will be affected by S-mediated attachment to 
9-O-Ac-Sias as well as by the close packing of HE molecules 
in the viral envelope. Serial dilutions of intact BCoV and OC43 vi- 
rions, adjusted for pNPA esterase activity (Table S2), were there- 
fore compared in solid-phase receptor-destruction assays with 
BSM (Figures 5C and 5D). Remarkably, the difference between 
the two viruses as measured by the onset of detectable receptor 
destruction was smaller than for their respective soluble HEs. 

The activity of virion-associated BCoV HE relative to that of 
soluble BCoV HE was reduced 8- to 16-fold. Conversely, the 
relative activity of virion-associated OC43 was increased at least 
in that destruction of receptors already became apparent at 
lower enzyme concentrations. (Figure 5C). To study whether 
these phenomena relate to the S protein and to S-mediated 
virion attachment to 9-O-Ac-Sias, we constructed recombinant 
murine coronaviruses (rMuCoVs) that express MuCoV S, which 
does not bind 9-O-Ac-Sias, in combination with either BCoV or 
OC43 HE. rMuCoV virions studded with OC43 HE lacked detect- 
able esterase activity toward BSM. In contrast, the esterase ac- 
tivity of rMuCoVs expressing BCoV HE was 6- to 8-fold higher 
than that of wild-type BCoV virions (Figure 5C). Apparently, 
B1CoV S proteins, by binding to 9-O-Ac-Sias negatively or posi- 
tively, affect HE esterase activity in BCoV and OC48 virions, 
respectively, depending on whether or not the HE lectin domain 
is functional. Most strikingly, BCoV and OC43 virions differed in 
the extent to which receptors were depleted. The soluble HEs of 
OC43 and BCoV, despite their difference in enzyme activity, de- 
stroyed all O-Ac-Sias as detectable by sp-LBA. BCoV virions de- 
stroyed 95% of the receptors after 24 hr incubation, which 
increased to 98% after 48 hr (Figure 5D). For OC43 virions, the 
effect was more pronounced such that after 24 hr only 85% of 
the receptors were destroyed. The remaining receptors were 
resistant even to continued incubation (Figure 5D). The data sug- 
gest that, in the context of the virion, the absence of a functional 
HE lectin domain results not only in reduced enzyme activity to- 
ward clustered receptors, but also prohibits cleavage of partic- 
ular receptor populations. As these “protected” receptor popu- 
lations are readily depleted by soluble OC43 HE (Figures 5B-5D), 
their resistance to cleavage by OC43 virions suggests that they 
are inaccessible to HEs when embedded in the viral envelope 
and in the presence of S. 

The data prompt a model in which the susceptibility of popu- 
lations of clustered receptors to cleavage depends on (1) their 
accessibility, as determined by the distance by which they 
extend from a fixed surface—the bottom of the ELISA well in 
our artificial system and, for example, the host cell membrane 
in vivo—and (2) the considerable difference in height between 
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Figure 5. Loss of HE Lectin-Mediated 
Receptor Binding Alters Sialate-9-O-Acety- 
lesterase Receptor, Destroying Activity 
toward Multivalent Substrates 

(A) Esterase activity of soluble recombinant 
HE*-Fc fusion proteins (SHE*) toward monovalent 
substrate pNPA. Enzymatic activity is shown as 
percentage of BCoV HE wild-type activity. 
WT, sHE* with wild-type HE ectodomains of 
BCoV-Mebus, OC43 USA/1967, or HCoV-HKU1 
as indicated; F2''A, BCoV-Mebus sHE* derivative 
with the lectin CBS inactivated through a 
FA substitution (Zeng et al., 2008); TREF, 
OC43 USA/1967 HE with the lectin CBS repaired 
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maxisorp plates, was incubated for 1 hr with 
2-fold serial dilutions of sSHE*. Receptor destruc- 
tion was assessed by sp-LBA with BCoV HE°-Fc 
at fixed concentration (2 ng/L). 

(C) Whole-virion receptor-destruction assays 
with 2-fold serial dilutions of purified viruses 
(starting at 15 mU). BCoV and OC43 sHE'*s were 
included for comparison. Receptor destruction 
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S and HE (20 and 8 nm, respectively), with the CBSs of S and 
those of HE esterase domains separated by a distance of up 
to 15 nm (Figure 5E). S-mediated virion attachment to a 
9-O-Ac-sialylated surface would bring the HE catalytic site in 
proximity of some clustered receptors (which would explain 
the stimulating effect of S on HE esterase activity in the case 
of OC43 virions, but not observed for rMuCoV-OC43 HE*) but 
at the same time would keep the virus-associated enzyme at 
“arm’s length” from other receptor populations, readily acces- 
sible to soluble HEs (which would explain the inhibiting effect 
of S on esterase activity in the case of BCoV virions and the 
apparent lack thereof in the case of rMuCoV-BCoV HE*) (Figures 
5C and 5D). Thus, a distinction arises between clustered recep- 
tor populations that come within immediate reach of the HE 
esterase and that therefore would be cleaved by virion-associ- 
ated HE as efficiently as by soluble HE and those that are kept 
out of reach of virion-associated HE as a result of S-mediated 


was measured after 1 hr or 4 hr incubation. Black 
dotted lines indicate 100% receptor destruction 
as determined with excess amounts of BCoV or 
OC43 sHE*. 

(D) Whole-virus-mediated receptor destruction 
over time with fixed amounts of BCoV-Mebus and 
O0C43 USA/1967 virions (15 mU). To exclude 
“exhaustion” (inactivation of OC43 esterase 
and/or changes in the physicochemical properties 
of the virions over time), we removed BCoV and 
OC43 at t = 24 and replaced them with equal 
amounts of freshly thawed aliquots of the virus 
stocks (black arrow). 

(E) Scaled side-by-side representations of the 
structures of coronavirus S (Walls et al., 2016) and 
HE proteins. Indicated are the estimated heights 
of S and HE and the approximate distance sepa- 
rating S receptor-binding sites and HE catalytic 
sites. (See also Figure S4 and S5 and Table S2.) 


receptor-binding and that are therefore cleaved at strongly 
reduced rates or even rendered non-cleavable unless HE is pro- 
vided with a functional lectin domain (Figure 6). By virions 
grasping on to multivalent sialoglycoconjugates via the HE re- 
ceptor-binding sites and “drawing them in” closer to the surface 
of the envelope, clustered receptors would yet become available 
to the HE esterase domain (Figure 6). Such a mechanism may be 
promoted by cooperativity among adjacent HE molecules and, 
membrane fluidity permitting, possibly even lead to redistribu- 
tion of surface projections in the viral envelope, resulting in 
displacement of spikes, local recruitment of HEs, and formation 
of a “receptor-destruction patch.” Saliently, some receptor pop- 
ulations may be completely resistant to destruction even by 
virion-associated HEs with a functional lectin domain, such as 
the 2% of detectable O-Ac-Sias remaining after 48 hr of incuba- 
tion with BCoV (Figure 5D). It is tempting to speculate that under 
natural conditions, the distinction between decoy and entry 
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Figure 6. Hypothetical Model for the Inter- 
action of HCoV and BCoV Virions with Multi- 
valent Glycoconjugates 

Schematically depicted are portions of HCoV 
(O0C43 or HKU1) and BCoV virions, with large 
spikes comprised of S (in gray) and smaller pro- 
trusions comprised of HE (yellow) extending from 
the viral membrane (blue). Functional HE lectin 
CBSs are indicated by black holes (only one 
shown per HE dimer for reasons of simplicity). Also 
shown schematically are membrane-anchored 
(bottom) and non-anchored (top) mucin-type 
glycoconjugates of bottle-brush filamentous 
appearance with clustered receptors (9-O-Ac- 
Sias, red dots) arranged in linear arrays and with 
absence of red dots indicating receptor-destruc- 
tion. The model, based on the size difference be- 
tween S and HE, visualizes how loss of HE lectin 
function might alter virion-associated receptor- 
destroying activity, reducing the specific activity of 
virions as well as the rate and selectivity of re- 
ceptor destruction. In virions with lectin-deficient 
HEs, clustered substrates will be largely kept at a 


distance from the HE esterase catalytic pocket as a result of S-glycoconjugate interaction. In contrast, HEs with intact lectin CBS may draw in portions of the 
glycoconjugates (or will draw the virion-associated HEs toward them), aided by cooperativity of binding between adjacent HEs within the viral envelope. Thus, 
clustered glycotopes become fixed within reach of the esterase catalytic sites and receptor destruction is promoted. 


receptors is made by default on basis of their accessibility to 
cleavage by virion-associated HE. 


CONCLUSION 


We showed that OC43—after its zoonotic introduction—has 
been under incessant selective pressure to adapt to the sialogly- 
come of the human host. HE-associated receptor binding was 
selected against, largely lost early on, and ultimately lost alto- 
gether through an accumulation of mutations in the lectin domain 
over a period of decades. In result, the balance between virion 
attachment and catalytically driven virion elution was reset, 
apparently to meet the specific requirements for optimal replica- 
tion in human airways and to allow the virus in this particular 
niche to distinguish between decoys marked for destruction 
and functional receptors that should be preserved for cell entry. 
This view is reinforced by our observation that HCoV-HKU1, 
another respiratory lineage A BCoV of humans, followed a 
convergent evolutionary path and also lost HE-mediated recep- 
tor binding. We offer that differences in the sialoglycomes of bo- 
vids and humans, such as in the local expression, structure, and 
density of 9-O-acetylated sialoglycoconjugates, may pose in- 
compatibilities during cross-species transmission of OC43 and 
BCoV and that the respective loss or preservation of HE-medi- 
ated receptor binding contributes to the host tropism of these 
viruses. This said, comparative sialoglycomics is still in its in- 
fancy, humans and bovids have not been studied for differences 
in O-Ac-sialoglycan expression in any detail, and the precise sia- 
loglycan-based constraints that selected for the particular traits 
of the animal and human BCoVs therefore remain to be identified. 
At any rate, our findings do reveal an as yet unappreciated 
aspect of lineage A BCoV adaptation to humans. Of broader 
relevance, they provide a general paradigm for dynamic, catal- 
ysis-driven virus-sialoglycan interactions in relation to host 
selectivity. This notion is supported by observations for influenza 
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A viruses, where neuraminidase (NA) activity toward multivalent 
substrates varies with differences in the length of the stalk 
domain and positively correlates with NA size (Castrucci and 
Kawaoka, 1993; Els et al., 1985), where NA catalytic activity is 
modulated by the absence or presence of a second Sia-binding 
site (Uhlendorff et al., 2009; Varghese et al., 1997) and where 
changes in NA stalk length and NA lectin affinity have been impli- 
cated in host specificity (Blumenkrantz et al., 2013; Castrucci 
and Kawaoka, 1993; Lai et al., 2012; Li et al., 2014; Matrosovich 
et al., 1999; Munier et al., 2010; Uhlendorff et al., 2009; Varghese 
et al., 1997). These observations can be readily interpreted in the 
context of our data and model proposed. 


EXPERIMENTAL PROCEDURES 


Materials and experimental procedures are detailed in the Supplemental 
Experimental Procedures. 


Expression and Purification of HEs 

Human codon-optimized sequences of the HEs of BCoV-Mebus, ECoV-NC99, 
PHEV-VW572, RbCoV-HKU14-1, CRCoV-240/05, HCoV-HKU1, and OC43 
USA/1967 (ATCC VR-759) were cloned in expression plasmid pCD5-T-Fc 
(Zeng et al., 2008). The resulting constructs encode chimeric proteins 
comprised of the HE ectodomain fused to the human IgG1 Fc domain, with 
the domains separated by a thrombin cleavage site. The fusion proteins 
were either expressed in an enzymatically active form (HE*) or rendered inac- 
tive through a catalytic Ser-to-Ala substitution (HE°). Plasmids coding for 
O0C43 NL/A/2005 and OC43 NL/B/2005 HE were constructed through site- 
directed mutagenesis of the OC43 USA/1967 HE expression vector using 
the Q5 kit (New England Biolabs). For lectin affinity and esterase activity as- 
says, HE-Fcs were produced by transient expression in HEK293T cells and pu- 
rified from cell culture supernatants by protein A affinity chromatography and 
low pH elution as described (Zeng et al., 2008). For crystallization, OC43 
NL/A/2005 HE*-Fce was expressed in HEK293 GnTI(—) cells and purified by 
protein A affinity chromatography, followed by on-the-beads thrombin cleav- 
age as described (Zeng et al., 2008). Beads were pelleted, and the HE ectodo- 
main in the supernatant was concentrated to ~15 mg/mL. HE genes from 
OC43 field strains were RT-PCR amplified with OC43 viral RNA, directly 


isolated from nasopharyngeal aspirates (van de Pol et al., 2006) as a template, 
and their sequences were used to construct codon-optimized versions. 


Purification of Virions for Whole-Virus Receptor-Destruction Assays 
Cell monolayers were infected at MOI 0.01. Virions were purified by ultracen- 
trifugation through 20% (w/v) sucrose cushions and resuspended in PBS. 
Virus preparations were analyzed for particle content by qPCR, plaque 
assay, quantitative latex bead ratio EM, and pNPA esterase activity assay 
(Table S2). 


Hemagglutination Assay 

HAA was performed with rat erythrocytes (Rattus norvegicus strain Wistar; 
50% suspension in PBS) and, in standard tests, with 2-fold serial dilutions of 
HE°-Fc proteins (starting at 25 ng/well) as described (Zeng et al., 2008). For 
high-sensitivity HAA, based on multivalency-driven high-avidity binding, 5 1g 
HE°-Fe was complexed to 2 x 10° protein A-coated 100 nm nanobeads 
(Chemicell) in 100 .zL PBS for 45 min at 4°C prior to 2-fold serial dilution. Hem- 
agglutination was for 2 hr at 4°C unless stated otherwise. 


Solid-Phase Lectin-Binding Assay and On-the-Plate O-Ac-Sia 
Depletion Assay 
sp-LBAs were performed with bovine submaxillary mucin (BSM) and 2-fold se- 
rial dilutions of HE°-Fc proteins as described (Langereis et al., 2012, 2015). 
Receptor-destroying esterase activities of soluble recombinant HEs (SHE+) 
and whole viruses were measured by on-the-plate O-Ac-Sia depletion assay, 
performed essentially as described (Langereis et al., 2015). BSM, coated on 
Maxisorp flat-bottom 96-well plates, was (mock-)treated with 2-fold serial di- 
lutions of SHE* in PBS (100 j:L/well) starting at 20 ng/L for 1 hr, or, for whole- 
virion receptor-destruction assays, with virions or sHE* starting at 15 mU 
esterase activity for up to 48 hr at 37°C. Depletion of O-Ac-Sia receptors 
was detected by sp-LBA with BCoV-LUN HE® (2 ng/L; 100 ptL/well). Esterase 
activity was plotted as the inverse of lectin-binding measurements, expressed 
in percentages. 
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